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We present electronic structure calculations, electrical resistivity data and the first specific heat 
measurements in the normal and superconducting states of quasi-one-dimensional A^MogSee (M 
= Tl, In, Rb). Rb2Mo6Seg undergoes a metal-insulator transition at ~ 170 K: electronic structure 
calculations indicate that this is likely to be driven by the formation of a dynamical charge density 
wave. However, T^MoeSee and In2Mo6Se6 remain metallic down to low temperature, with super- 
conducting transitions at T c = 4.2 K and 2.85 K respectively. The absence of any metal-insulator 
transition in these materials is due to a larger in-plane bandwidth, leading to increased inter-chain 
hopping which suppresses the density wave instability. Electronic heat capacity data for the su- 
perconducting compounds reveal an exceptionally low density of states De f — 0.055 states eV _1 
atom" 1 , with BCS fits showing 2A/k B T c > 5 for Tl 2 Mo 6 Se 6 and 3.5 for In 2 Mo 6 Se 6 . Modelling 
the lattice specific heat with a set of Einstein modes, we obtain the approximate phonon density 
of states F(lj). Deconvolving the resistivity for the two superconductors then yields their electron- 
phonon transport coupling function oq T F{ui). In T^MoeSeg and In2MoeSe6, F{uj) is dominated by 
an optical "guest ion" mode at ~ 5 meV and a set of acoustic modes from ~ 10-30 meV. Rb2MoeSe6 
exhibits a similar spectrum; however, the optical phonon has a lower intensity and is shifted to ~ 
8 meV. Electrons in Tl2Mo6Se6 couple strongly to both sets of modes, whereas In2MoeSe6 only 
displays significant coupling in the 10-18 meV range. Although pairing is clearly not mediated by 
the guest ion phonon, we believe it has a beneficial effect on superconductivity in T^MogSee, given 
its extraordinarily large coupling strength and higher T c compared to In2MoeSe6. 

PACS numbers: 71.30.+h, 74.25.-q, 74.70.Dd 



I. INTRODUCTION 

The M 2 Mo 6 Se 6 (M = Tl, In, Rb, Li, Na, K, Cs) sys- 
tem was first discovered by Potel et a/P and is closely 
related to the well-known quasi-three-dimensional (quasi- 
3D) Chevrel Phase compounds. Rather than comprising 
individual "zero-dimensional" MoqXs (X = S, Se, Te) oc- 
tahedral clusters coupled by a metallic cation, these ma- 
terials are composed of quasi- ID (MogSee)^ chains ori- 
ented along the z axis, weakly coupled by M ions. Only 
T^MogSeg and In 2 Mo6Se6 are superconducting, with T c 
= 3 - 6.5 K (varying between samples) and ~ 2.9 K re- 
spectively. In contrast, Rb 2 Mo6Se6 undergoes a broad 
metal-insulator transition between 100 K and 200 K. 2 
Little data currently exists in the literature for the re- 
maining members of the family, although it is known 
that they exhibit similar metal-insulator transitions and 
do not become superconducting under ambient pressure 
at low temperature!^ 

Reduced dimensionality and its effect on supercon- 



ductivity remains one of the central issues in contem- 
porary condensed matter physics research. Since the 
late 1970s, numerous unconventional superconductors 
displaying highly anisotropic properties in both the nor- 
mal and superconducting states have been discovered. 
Among these, notable examples include the quasi-2D 
high-temperature cuprate superconductors (HTS) and 
the quasi-lD organic Bechgaard salts. However, the 
most strongly ID superconductors synthesized to date 
have attracted remarkably little attention over the years. 
Tl 2 Mo 6 Se 6 and In 2 Mo 6 Se6 boast anisotropy ratios e = 

H&IHk ^ 12 -° and 17 - 2 respectively, 1 ^ significantly 
greater than e « 8.5 in (TMTSF) 2 C10 4 P Furthermore, 
these materials do not possess any intrinsic magnetism, 
thus rendering them an ideal uncomplicated system for 
the study of low-dimensional superconductivity. 

In comparison with the HTS, the number of publi- 
cations existing for M 2 Mo6Se6 is around three orders 
of magnitude sma ller. Early work concentrated on the 
electrical transport!^ and magnetic properties,^ imme- 



diately revealing large anisotropies in both the normal- 
state resistivity and superconducting coherence length 
for Tl 2 Mo 6 Se 6 . Two distinct classes of Tl 2 Mo 6 Se 6 were 
identified 6 -' by the behaviour of their longitudinal resistiv- 
ity: A-tyj>e samples with conventional metallic behaviour 
down to low temperature or 5-type samples displaying 
a broad minimum for T < ~ 80 K followed by an up- 
turn reminiscent of charge density wave (CDW) forma- 
tion. More recent measurements^ have shown that the 
coherence length perpendicular to the chain axis £j_ is 
at most 75 A, a value not significantly larger than that 
found in some HTS. In an important parallel with organic 
quasi-lD superconductors,^ the upper critical field per- 
pendicular to the chain axis z, H^, does not saturate 
down to 50 mK. 12 Hall effect measurements by the same 
authors display a regime crossover at T ~ 80 K which 
they attribute to the onset of a CDW or spin density 
wave (SDW). However, there is no support for the forma- 
tion of a CDW in normal-state resistivity curves for the 
A-type samples, which nonetheless display a Hall effect 
crossover. Furthermore, the weak temperature-invariant 
diamagnetism in Tl 2 MogSeg revealed by ac susceptibil- 
ity measurements^ does not encourage a SDW interpre- 
tation. 

The discovery that the application of uniaxial stress 
along the z axis in Tl 2 MogSe6 suppresses superconduc- 
tivity and induces a metal-insulator transition increased 
the evidence for this system being close to a CDW insta- 
bility, in particular due to the non-linear I— V curves and 
broadband noise characteristic of density wave motion 
observed.^ Conversely, hydrostatic pressure increases 
the conductivity of Tl 2 MoeSe6 in the normal state but 
still suppresses superconductivity. 14 

At first glance, the quasi-one-dimensional nature of the 
il^MogSee family renders them strong candidates to un- 
dergo a Peierls transition, so it was initially a mystery 
as to why Tl 2 MogSe6 and In 2 MogSe6 remained metallic 
down to low temperature. Early attempts to resolv e this 
question focussed on the calculated band structur e! 15 * 16 ! 
which displayed three contributions to the Fermi sur- 
face: a broad singly-degenerate quasi-lD Mo d "helix" 
band and two doubly-degenerate 3D "octahedral" elec- 
tron pockets at the zone boundary. The occupancy of 
these 3D pockets was believed to stabilise the structure 
against a Peierls transition and the authors of all ensuing 
publications attempted to interpret their results within 
this multi-band framework. However, band structure cal- 
culation techniques have significantly advanced since the 
non-self-consistent approaches of the early 1980s. It is 
therefore instructive to recalculate the band structure 
of M 2 MoeSe6 using a fully self-consistent method and 
compare our results with the existing calculations for 
Tl 2 Mo 6 Se 6 . 

The nature of the superconducting state in Tl 2 Mo6Se6 
and In 2 MoeSe6 also remains very unclear, particularly 
since recent scanning tunnelling microscopy (STM) ex- 
periments^ on the related quasi-3D Chevrel Phase su- 
perconductor PbMogSs have provided strong evidence 



for a highly anisotropic or noded gap function. In- 
elastic neutron scattering measurements of the phonon 
density of states (PDoS)^ revealed a strong low-energy 
Einstein-like optical mode attributed to vibrations of the 
M "guest" ions between the MogSeg chains, as well as 
higher energy intra-chain modes similar to those seen in 
the 3D Chevrel phases P^ This finding is supported by 
early normal-state specific heat data 2 P which unfortu- 
nately lacked sufficient resolution to provide any infor- 
mation on the superconducting state. However, no stud- 
ies of the electron-phonon coupling were performed; nor 
have any tunnelling experiments been carried out on the 

Af 2 MogSe6 system. 

Motivated by recent discoveries in boridc42M2J anc [ ft. 
pyrochlore 2 ^ systems where superconductivity is medi- 
ated by a low-energy rattling phonon, we therefore de- 
cided to measure the PDoS and electron-phonon coupling 
for Tl 2 MogSe6, In 2 Mo 6 Se6 and Rb 2 MoeSe6 by decon- 
volving normal-state specific heat and resistivity data. 
These measurements are analysed in parallel with our 
specific heat and resistivity data below T c for M = Tl,In 
and in the insulating phase for M = Rb. In addition, 
we have performed a complete theoretical analysis of the 
electronic structure and parameters governing the metal- 
insulator transition in Af 2 MoeSe6. By combining our ex- 
perimental data with the trends indicated by our new 
band structures, we hope to remove some of the con- 
fusion surrounding the mechanism for superconductivity 
and the effects of low-dimensionality in this fascinating 
system on the borderline between superconducting and 
insulating instabilities. 



II. THEORY 

A. Crystal Structure 

The crystal structure of the compounds with chemical 
formula M 2 MogSe6 is shown in Fig. [l] viewed both par- 
allel and perpendicular to the z axis. Mo and Se atoms 
form quasi-lD (Mo6See)oo chains oriented along the z 
axis, separated by M ions in a zig-zag formation. The 
chains consist of inner Mo and outer Se triangles, stacked 
with a c/2 separation along the z axis and rotated 180° 
with respect to each other. The axes of the Mo-Se trian- 
gles are aligned with each other, but rotated by 10° from 
the lattice vectors in the x — y plane. We may consider 
the chains to be a linear condensation of MoeSs clusters 
via face-sharing of the Mo6 octahedra. These clusters 
are the building blocks of the related quasi-3D Chevrel 
phases. 

The conventional unit cell is hexagonal (space group 
P63/W), and contains two formula units (f.u.)} M atoms 
occupy 2d positions while Mo and Se atoms occupy 6h 
positions. The hexagonal lattice parameters (determined 
by X-ray diffraction) arc a H = 8.854, 8.934, 9.257 A and 
c H = 4.493, 4.494, and 4.487 A respectively for In, Tl 
and Rb-based crystals. It can immediately be seen that 
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FIG. 1: Crystal structure of A^MogSes compounds in the 
001 plane (above) and side view (below) of the (MogSee),^ 
chains. Darker (lighter) symbols indicate atoms sitting on 
even (odd) planes respectively. The local coordinate system 
used to plot the partial Mo, Se characters in Fig. [2] is also 
shown. 
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In 
Tl 
Rb 



%Mo 



VMo 



XSe 



VSe 



0.189 
0.187 
0.181 



0.156 
0.154 
0.149 



0.068 
0.067 
0.064 



0.369 
0.366 
0.355 



TABLE I: Optimized internal coordinates (Wyckoff positions) 
for Mo and Se in the three M2Mo6Se6 compounds considered 
in this work (M=Tl,In,Rb). 



the inter-chain distance correlates with the atomic radius 
of the M atom, whereas Ch (and hence the intra-chain 
atomic separation) remains roughly constant regardless 
of M . Furthermore, the intra-chain atomic separations 
are much smaller than the inter-chain distances. For ex- 
ample, in Tl the shortest Mo-Mo distance (for two Mo 
atoms in the same triangle) is 2.66 A; the shortest Mo-Se 
distance is 2.695 A and the Se-Se separation is 3.767 A. 
In contrast, the Mo-Mo inter-chain separation is 6.34 A. 
The crystal structure thus provides an immediate in- 
dication of the strong anisotropy present in this family 
of materials. In the following, we will show that a large 



anisotropy is also found in the electronic structure and 
discuss its consequences for the physical properties of 
Af 2 Mo 6 Se 6 compounds. 



B. Electronic Structure 

We have performed ab — initio Density Functional 
Theory (DFT) calculations of the electronic properties 
of M 2 Mo 6 Se 6 , for M =T1, In, Rb, employing the the 
full-pot ential Linear Augmented Plane Wave (LAPW) 
method!- 24 * 25 ^ For all systems, we used the experimen- 
tal lattice parameters and optimized the internal coor- 
dinates: these optimized values are given in table IT] 

For T^Mc^Se^ where the experimental Wyckoff posi- 
tions are knownj 2 ^ the optimized coordinates agree with 
the experimental ones to better than 1 %; for R^MogScg 
and Rb2MogSe6 we could not compare with experimen- 
tal data. Structural optimization yielded intra-chain dis- 
tances which do not depend (within computational accu- 
racy) on the nature of the M atom, whereas interchain 
distances increase by 5 % going from Tl,In to Rb due to 
the larger in-plane lattice constant. 
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FIG. 2: Band structure and partial DoSs of TbMoeSee; en- 
ergies are in eV and measured w.r.t. the Fermi level, DoS 
are in states/ (eV cell). Colors indicate different I charac- 
ters (green=s, orange=p, red=d). Shaded areas, delimited by 
dashed lines, refer to the partial characters (Tl p z , Mo d xz , 
Se p x ) which provide the largest contribution at E_p. For a 
definition of the axes and atoms, see Fig. [I] An enlargement 
of the DoS around the Fermi level is also shown; the numbers 
indicate the enlargement factor. 

Since the major features of the electronic structure 
of M2MoeSe6 are the same for M=T1, In, Rb, we will 
first present the electronic structure of Tl 2 Mo 6 Seg and 
then discuss the differences with In and Rb. Our calcu- 
lated electronic structure reproduces the m ain fe atures 
of the earlier non-self-consistent calculations] 15 * 16 * except 
for one crucial difference in the position of the Fermi 
level. The electronic bands and partial densities of states 



(DoS) are shown in Fig. [51 with all energies measured 
relative to the Fermi level. Se s and Tl d states form 
weakly-dispersing bands, located 15 eV and 12 eV below 
the Fermi level respectively. The Tl 6s states lie approx- 
imately 7 eV below the Fermi level; Tl is therefore in 
a nominal +1 valence state, effectively behaving like an 
alkali metal. The bands lying between -5 and +5 eV are 
mainly of Mo d and Se p character, although there is a 
significant hybridization with Tl p states as can be seen 
in Fig. [2] The Mo sp bands are displaced to higher ener- 
gies by covalent hybridisation with Se p bands and hence 
lie above +5 eV. 

The 18 Se p bands are centered around -4 eV and over- 
lap with the lower portion of the 30 Mo d bands, which 
extend over ±5 eV around the Fermi level. These are 
divided into 12 bonding and 16 antibonding states, sep- 
arated by a pseudogap ~ 1 eV wide around the Fermi 
level. Two highly ID bands cross the gap: they de- 
rive from the zone-folding of a single helix band. The 
Fermi level cuts the band structure at k z =ir/c, exactly 
at half filling. The corresponding Fermi Surface is shown 
in Fig. [3j In contrast with earlier calculations, we do not 
find the additional small 3D pockets at the Fermi level 
which were previously alleged to be responsible for the 
stability of the chains against any Peierls distortion and 
ensuing density wave transition.^ The octahedron band 
(of Mo d 3z 2 — r 2 and x 2 — y 2 character) responsible for 
these pockets in fact lies ~ 0.1 eV above Ep. Fig. [4] 
shows that In 2 Mo6Se6 and F^MogSeg have a very sim- 
ilar band structure to T^MogSeg, i.e. there are no other 
bands at Ep except for the helix band. Small variations 
in the dispersion of this band must therefore control the 
stability of the Af 2 Mo6Se 6 family. 

Let us now consider the relative anisotropy of the 
M2M06SC6 series, in terms of details of the dispersion 
of the helix band, using a simplified tight-binding model 
with large out-of-plane (W) and small in-plane (w) band- 
widths. In the three rightmost panels of Fig. [2] we high- 
light the partial orbital characters which give the largest 
contribution to the helix band: the dominant contribu- 
tion is from Mo d xz states. (The labels of the orbitals 
refer to the local coordinate systems centered on the Mo 
atoms, shown in Fig. fl]). 

In fact, the Bloch states that form this band are built 
from the in-phase linear combination of the three equiva- 
lent Mo xz orbitals which sit on each M03 triangle. Due 
to the presence of a two-fold screw axis along the cen- 
ter of the chain, they form right and left-handed helices 
which wind up along the chain axis. (For a definition of 
the axes, see Fig. [I]). The out-of-plane bandwidth W is 
thus given by the hopping of one Mo d xz to another Mo 
d xz orbital on the next plane; the in-plane bandwidth is 
given by the hopping between the chains, mediated by Se 
p x and Tl(In) p z or Rb s orbitals (illustrated by yellow 
and green lines respectively in Fig. [2]) . 

Due to symmetry, hopping via Tl(In) p z orbitals is 
more effective than via Rb s orbitals. In Fig. [4] we show 
a close-up of the band structures of A^MoeSee, decorated 
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TABLE II: Parameters governing the metal-insulator transi- 
tion in A^MosSee [M = Tl,In,Rb) from local density approxi- 
mation (LDA) calculations. W and w are the out-of-plane and 
in-plane band widths of the helix band in eV: they determine 
the dimensionless critical coupling constant A c (see Eq. [2b. 
Wo,» are the bare frequencies of the two Peierls modes in meV. 
A(a>o,0 are the corresponding dimensionless bare electron- 
phonon coupling constants (see Eq. [lj; 2 So is the value of 
the Peierls gap (Eq. E| in eV. The renormalized phonon fre- 
quencies u)i (in meV) are obtained from Eqs. [3] and [5] for Tl, 
In and Rb respectively. 
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FIG. 3: Fermi surface of Tl2Mo6Se6, in a vertical plane cut- 
ting through the center of the hexagonal Brillouin zone, shown 
in Fig. [I| 



with the partial characters associated with the M atoms. 
Here, we observe that the contribution of Tl and In p z 
states to the helix band is much stronger than that of the 
Rb s states. This agrees well with the significant In p or- 
bital contribution to the conduction band in In 2 Mo 6 Sc 6 
suggested by nuclear magnetic resonance spectroscopy.^ 
Correspondingly, the in-plane bandwidth is reduced by a 
factor of ~ 10 in Rb2MogSe6. From a fit of these three 
band structures, we extract the values of W and w re- 
ported in table [IT] 




FIG. 4: Band structures of M 2 Mo 6 Se 6 , for M=Tl,In and 
Rb decorated with partial M characters. Energies, measured 
w.r.t. the Fermi level, are in eV. 



C. Mechanisms for Metal-Insulator Transitions 

The first possibility for the metal-insulator transition 
is a CDW instability, due to the interaction of the he- 
lix band with phonons (i.e. a Peierls transition). In 
Af 2 Mo 6 Se6, there are two Peierls-like modes which could 
lead to dimerization of the chains and open a gap at Ep : 
the eigenvectors are such that the Mo triangular units 
on subsequent planes move out-of-phase with respect to 
each other and the surrounding Se triangles can then be 
displaced either in or out of phase. 

Examining these two modes, we use local density ap- 
proximations (LDA) to determine two bare frequencies 
wo,i and u>o,2 which we have found to be independent of 
the nature of the M guest ion. Our calculated values are 
Wo,i = 22 meV and aj . 2 = 27 mcV. These frequencies in- 
clude the response of the whole electronic band structure 
to the phonon perturbations, with the exception of the 
helix band, which we calculate analytically below. This 
procedure is more accurate than a direct LDA calcula- 
tion of the phonon frequency, which becomes extremely 
inaccurate as the in-plane bandwidth w of the helix band 
tends to zero. 

The effect of both modes on the band structure is 
the same: a frozen displacement eu of the atoms along 
the renormalized eigenvector opens a gap 25 = 2Du at 
the Fermi wavevector k z — ir/c. This defines the de- 
formation potential D which (together with the bare 
phonon frequencies ujQ,i and the DoS at the Fermi level 
N(0) = 1/W) determines the bare electron phonon cou- 
pling constants A(wo,i) through the relation: 



A(wo,i 



D 2 



wMui; 



(i) 



Due to the finite in-plane bandwidth of the helix band, a 
Peierls transition can take place only if this bare electron- 
phonon coupling constant A(wo ») exceeds a critical value 
A,: 



(2) 



If A(wo,») < A c , there is no Peierls transition, but the 
bare phonon frequency is renormalized: 



Wi =w ,h 1 - 



A(uJo,i) 



(3) 



If A(u>o,i) > A c , there is a Peierls transition, with a gap 
2(5 = 2Du : 



2S = 2Wexp(- 



I 



2A(w ,i) 



)• 



(4) 



The frequency of the i phonon in the off-centre posi- 
tion is then given by: 



: w ,iy 2A(w ,i)- 



(5) 



2 (hx(W/w) + 1) ' 



Table |TT] lists the relevant LDA parameters for 
M=Tl,In,Rb. It may readily be seen that all ilf 2 Mo6Se6 
compounds are close to the CDW instability, since the 
values of the critical interaction parameters A c are very 
low compared to the bare electron-phonon coupling pa- 
rameters Ao.i. 

For T^MogSeg and In 2 MoeSe6, the bare electron- 
phonon coupling constants for both the lower and 
the higher Peierls modes are slightly below A c ; for 
Rb 2 Mo 6 Se 6 A(u>o,i) is well above the critical value, thus 
indicating that a Peierls transition is likely to occur with 
a gap 2(5 = 0.34 eV. However, due to the small value of 
A(wo,i)) the off-center minimum of the phonon potential 
is very shallow, implying that the transition is of dynam- 
ical rather than static character. This means that the 
density wave formation may not be characterised by a 
static structural modulation, as is the case for "classi- 
cal" CDW systems such as NbSe2. 

An alternative explanation for the observed metal in- 
sulator transition in M 2 Mo6Se6 compounds could be a 
spin density wave (SDW) instability. A staggered an- 
tifcrromagnetic order of M03 units is the spin analogue 
of the Peierls distortion. The spin density wave (SDW) 
transition is regulated by the magnetic coupling constant 
A/=||^f, where the Stoner parameter / can be esti- 
mated from the atomic value for Mo: 0.6 eV. A/ is one 
order of magnitude smaller than A c in Tl, In and Rb; we 
therefore do not find a stable antiferromagnetic (AFM) 
solution by LDA for any of the M 2 Mo 6 Se 6 compounds. 

However, AFM ordering is also favored by the Coulomb 
repulsion U, which is much larger than Hund's cou- 
pling /. Constrained LDA calculations give U—5.1 eV 
and J=0.62 eV per Mo atom in Tl 2 Mo 6 Se 6 . With 
these values, we estimate a magnetic coupling constant 
A;7=0.12 for Tl,In and 0.11 for Rb, which implies that all 
MaMogSeg are on the verge of a SDW transition due to 
strong electronic correlations. 

For both charge (Peierls) and spin (AFM) density 
wave transitions, the stronger tendency of Rb 2 Mo 6 Se 6 
towards insulating behaviour compared with Tl 2 MogSeg 



and Iri2Mo6Se6 is due to the smaller value of the criti- 
cal coupling constant A c . This is mainly caused by the 
reduced in-plane bandwidth of the helix band, which we 
recall is one order of magnitude smaller for M = Rb 
than for M = In, Tl. Further details of our calcula- 
tions and the analytical model employed will be discussed 
elsewherePS We will now turn to our experimental data, 
which confirms our experimental calculations and pro- 
vides further insight towards the metal-insulator transi- 
tion. 



III. EXPERIMENT 

A. Sample Preparation and Experimental Details 

Needle-like crystals of dimensions approximately 
4 mm x 300 /^mxlOO /xm and mass ~ 800 /^g were syn- 
thesised by different methods depending on the thermal 
stability of the compounds: TL^MoeSee and Rb 2 MogSe6 
were prepared in sealed molybdenum crucibles at 1700 °C 
and 1500 °C respectively, whereas In 2 Mo 6 Se 6 was pre- 
pared in an evacuated sealed silica tube at 1100 °C. Their 
crystalline structures were verified by the mono-crystal 
diffraction method using a KAPPA CCD NONIUS and 
exhibit a slight cation deficiency with occupancy factors 
0.95, 0.94 and 0.93 for In, Tl and Rb based compounds re- 
spectively. Larger polycrystalline samples of Tl 2 MogSe6 
and In 2 MogSe6 of mass 10-75 mg were also produced by 
the same methods. 

The samples were initially characterized by AC Suscep- 
tibility, measured using a Quantum Design™ Physical 
Property Measurement System (QD PPMS). As shown 
in Fig. [5j Tl 2 Mo 6 Se 6 and in 2 Mo 6 Se 6 both exhibit su- 
perconducting transitions at T c =4.2 K and 2.85 K, with 
AT C =0.7 K and 0.5 K respectively. Rb 2 MogSe 6 displays 
a constant weak diamagnetic signal down to the lowest 
temperature measured (1.7K) with no sign of supercon- 
ductivity being observed. No significant increase in tran- 
sition width was observed for the polycrystals compared 
to the single crystals, indicating a high sample quality. 
We stress that all members of the M 2 Mo6Se6 family re- 
main stable under atmospheric conditions and none of 
our samples has exhibited any ageing effects. 

AC resistivity was measured using the same QD PPMS 
with a Helium-3 insert from 0.35-300 K. Four gold con- 
tacts of thickness ~5 nm were sputtered onto single crys- 
tals of each compound and 50 /mi gold wires attached 
using silver epoxy glue. This method yielded contact re- 
sistances ~1 ft. Short (1 s) pulses of a small AC current 
(0.02 mA, 470 Hz) were used to minimise any heating 
effects in the sample at low temperature. 

Specific heat was initially also measured in the QD 
PPMS using a standard relaxation technique from 0.35- 
300 K. The largest homogeneous polycrystalline samples 
available for each compound were mounted using a mea- 
sured quantity of Wakefield™ grease (whose contribu- 
tion to the heat capacity was later subtracted). How- 
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FIG. 5: AC susceptibility of Tl 2 Mo 6 Se 6 , In 2 Mo 6 Se 6 and 
Rb2Mo6Se6 from 1.8-300 K. Inset: zoom onto superconduct- 
ing transitions, with T c =4.2 K and 2.85 K in Tl 2 MogSe6 and 
I^MogSeg respectively. As is the case for all graphs in this 
work, only a small fraction (typically 10-20 %) of the data- 
points measured are explicitly marked for clarity. 



ever, due to the extremely small density of states at the 
Fermi level the QD PPMS was unable to detect the su- 
perconducting transition in Tl 2 MogSe6 and In 2 MogSeg. 
High-sensitivity relaxations from 1.3-10 K were therefore 
carried out on single crystals in our dedicated specific 
heat laboratory, enabling us to study the superconduct- 
ing transition with a magnetic field both perpendicular 
and parallel to the z axis. 



B. Characteristics of the Superconducting State in 
Tl 2 Mo 6 Se 6 and In 2 Mo 6 Se 6 

As previously reported,^ resistive transitions into the 
superconducting state of Tl 2 MogSe6 exhibit an anoma- 
lous broadening under an applied magnetic field H, sim- 
ilar to that seen in the HTS. We have performed similar 
measurements on In 2 Mo6Se6 and show these in Fig. |6j 
together with the results in Tl 2 MogSe6 for comparison. 
It can immediately be seen that In 2 MogSe6 also displays 
a broadening of the transition with increasing field, al- 
though the effect is less spectacular than in Tl 2 Mo 6 Se 6 . 
This is highlighted in Fig. 6jf), where we have plotted the 
normalised transition width AT C (H) / AT c (0) as a func- 
tion of normalised perpendicular magnetic field H/H c % 
for both Tl 2 Mo6Se6 and In 2 MogSe6. Both materials dis- 
play a linear behaviour in AT C (H) / AT c (0) as the field 
increases, with the gradient for Tl 2 MogSe6 ~ 30% larger 
than that of In 2 MogSe6. 

The width of a superconducting transition is gov- 
erned by two parameters: the intrinsic homogeneity of 
the superconductor and the narrow thermal fluctuation- 
dominated critical region which surrounds any phase 
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FIG. 6: (a)-(e): Resistive transitions in TbjMoeSee and 
In2Mo6See with field both parallel and perpendicular to the 
z axis. Two samples of T^MoeSeg are shown: Sample (1) 
displaying a double transition indicating an inhomogeneous 
Tl content and Sample (2) with a single broad transition, (f): 
Normalised transition widths AT C (H) / AT c (0) as a function 
of normalised magnetic field H/H C 2. For a transition in mag- 
netic field H, AT C is defined as the temperature difference 
between resistivities of 5% and 95% of the saturated normal- 
state value at 7 K. 



transition. Inhomogeneity contributions are field- 
independent, but for 3D fluctuations the critical region, 
whose width is given by the Ginzburg parameter G^d = 
{k b T c /V2£i£\\H c (0) 2 ) 2 multiplied by T c , is weakly field- 
dependent, since G 3D (H) « G\^{Hj 'H c2 ) 2 / 3 . This 
clearly cannot explain our data, since we observe AT C 
to increase linearly as H increases. However, Mishonov 
et auF®* derived a quasi-lD Ginzburg parameter for a su- 
perconducting nanowire: 



G\d — 



kj 



8V^AC£(0)S 



(6) 



where AC is the specific heat jump at T c , S is the 
cross-sectional area of the nanowire and we only con- 
sider longitudinal fluctuations so £(0) = £||(0). Modelling 
A^MoqScq as a weakly-coupled assembly of supercon- 
ducting filaments, each with the radius of a single MogSeg 
chain, we calculate critical region widths G\dT c = 1.5 K 
and 2.0 K for Tl 2 Mo6Se6 and In 2 Mo6Se 6 respectively. 

Similar trends are also seen in the electronic contri- 
bution to the heat capacity (Fig. [7]), with Tl 2 Mo 6 Se 6 
and In2MogSe6 behaving very differently in a magnetic 
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FIG. 7: (a)-(c): Superconducting transitions in T^MosSee 
and In2Mo6Se6 seen by specific heat with field both parallel 
and perpendicular to the z axis. (d),(e): BCS fits of the zero- 
field transition in TFjMosSee and I^MoeSeg. 



field. The onset temperature T ons of the specific heat 
jump in In 2 Mo6Se6 is rapidly and uniformly displaced 
to lower temperature with increasing H, whereas regard- 
less of the applied field strength or orientation, T ons does 
not drop below ~ 4K in T^MogSeg. It should be noted 
that In2MoeSe6 shows an abnormally large electronic 
contribution to the specific heat 7 in proportion to its 
specific heat jump, suggesting that only around 50% of 
the electrons at the Fermi level become superconducting. 
This scenario could be explained by a slight variation in 
the sample stoichiometry leading to a CDW coexistent 
with superconducting regions in the same crystal, hence 
also supporting our observation that In 2 Mo6Se6 is more 
anisotropic than T^MogSeg. Correcting for this anomaly 
in In 2 Mo6Se6, we estimate a Sommerfeld constant 7 « 
0.13 mJgat _1 K- 2 for both Tl 2 Mo 6 Se 6 and In 2 Mo 6 Se 6 , 
corresponding to a dressed density of states at the Fermi 
level De f — 0.055 states eV _1 atom -1 . 

We have attempted to fit the low-temperature spe- 
cific heat in Tl 2 Mo6Se6 and In 2 Mo 6 Se 6 using a standard 
BCS s-wave single-band a-model, as shown in Fig.^d). 
As discussed in the previous section, our band struc- 
ture calculations indicate that only a single band (the 
ID Mo d helix) crosses the Fermi level, thus eliminat- 
ing any possibility of multi-band superconductivity in 
A/ 2 MogSe6. Although the jump we measure at T c is ther- 
mally broadened in each case, this yields an excellent fit 
for In 2 MoeSe6 with a gap value of 0.4 meV correspond- 
ing to 2Ao/fcsT c = 3.5, just below the standard weak- 
coupling BCS value 3.52. An unusual hump on the back 
of the peak in Tl 2 Mo 6 Seg renders the fit more difficult 
in this compound, but by using entropy conservation it 
is still possible to estimate a gap Ao > 1.1 meV and 
2Ao/ksT c > 5, placing it in the extreme strong-coupling 
regime. Although the quality of our fit in In 2 Mo 6 Se 6 
does not appear to favour the presence of low-lying ex- 
citations, we were not able to accurately measure to suf- 
ficiently low temperatures in Tl 2 MoeSe6 to conclusively 
rule out the existence of <i-wave superconductivity in this 
material. However, it is clear that Tl 2 MoeSe6 has a sig- 
nificantly less conventional superconducting ground state 
than In 2 Mo 6 Se 6 , an astonishing difference considering 
the close similarity between the compounds. 

Completing the analysis of our data using anisotropic 
Ginzburg-Landau theory,^ we summarise the supercon- 
ducting parameters of both Tl 2 Mo 6 Se6 and In 2 Mo6Se 6 
in Table |IH| The large calculated values for k highlight 
both the extreme type II nature of these superconductors 
and their enormous anisotropy. H^ 2 and £" should be re- 
garded as minima, due to the high sensitivity of these 
materials to the field orientationP We estimate a sample 
alignment better than ± 2 ° with the field; however non- 
parallel crystalline intergrowths may exist within a single 
needle-like sample which would reduce our measured H^ . 

The measured H\ 2 — 4.35 T in In 2 Mo 6 Se6 approaches the 
Clogston limit Hp = 4.9 T and it would be instructive 
to re- measure the resistivity of a small single crystal in a 
high-accuracy goniometer in order to verify the possibil- 



TABLE III: Measured and derived anisotropic superconduct- 
ing parameters in T^MogSee and I^MoeSeg 

Tl 2 Mo 6 Se 6 In 2 Mo 6 Se 6 

|| -L || -L 

Measured 

T c 4.2 K 2.85 K 

H c2 (0) 5.9 T 0.47 T 4.35 T 0.25 T 

H c (0) 0.0207 T 0.0119 T 

Derived 

e 12.6 17.2 

£(0) 940 A 75 A 1500 A 87 A 

A(0) 0.12 /mi 1.5 /im 0.13 /ton 2.2 /im 

K 202 1.3 260 0.87 

G 3D 3.3 10~ 6 3.0 10~ 6 

Gin 0.36 0.69 

7„ 0.13 mJ K -2 gat -1 0.13 mj K -2 gat -1 

BCS gap A > 0.9 meV 0.4 meV 

2A /k B T c _ > 5 3.5 

H P = Ao/v^Ms 11 T 4.9T 



ity of Pauli-limited superconductivity occurring in this 
compound. 



Normal-State Specific Heat and Phonon 
Densities of States 



Two features are immediately apparent in a plot of 
Ctot/T vs. T: a strong peak at T^ 80 K and a shoulder 
at T~ 20 K, indicating two dominant phonon energies. 
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FIG. 8: Total specific heat divided by temperature C ' /T for 
Tl2Mo6Se6 (circles), I^MogSee (squares) and Rt^MosSee 
(diamonds). Note that the increased noise level above 100 K 
in Rb2Mo6See can largely be attributed to a low sample mass. 
Inset: Close-up of C(T) in Rb2Mo6See, showing a small peak 
between 170 K and 185 K. 

The high-temperature specific data measured on larger 
polycrystalline samples are sufficiently noise-free to per- 



mit an inversion of the phononic contribution C p h(T) to 
the total heat capacity, hence obtaining the PDoS F(co). 
We stress that this method does not provide a detailed 
PDoS map of the type revealed by neutron scattering, but 
rather produces a smoothed phonon distribution func- 
tion F(u>). The specific heat and low-temperature fea- 
tures of the PDoS are accurately reproduced by F(u>), 
but it does not give a high-precision representation of 
the PDoS at high temperature. We model F(ui) as a se- 
ries of logarithmically-spaced Einstein modes with fixed 
energies and adjustable weights: 



then given by 
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FIG. 9: (a) Lattice specific heat divided by temperature 
G'/T — 7 for Li2Mo6Se6 (triangles) fitted (solid line) by an as- 
sembly of Einstein modes (dashed lines), (b) Lattice specific 
heat normalised by temperature cubed C — jT/T (triangles), 
highlighting the quality of our low temperature fit (solid and 
dashed lines) and the dominant contribution from the optical 
mode generated by the In ion. 



F(u) = J2 F kS{u-u k ) 



(7) 



C(T) = 3RY,Fk 



P^ ci-^k 



( e x k _ 1)2 



(8) 



where Xk = w^/T and ujk+i/^k = 1-75 to limit the 
number of modes and provide a stable solution. A least- 
squares fit of the measured specific heat for each com- 
pound was performed and the decomposition into Ein- 
stein modes shown for I^MogSeg as an example in Fig. [9] 
The fitting technique for T^MoeSee and PJ^MogSeg is 
identical and the results of a similarly high quality, accu- 
rately reproducing our experimental data over the entire 
temperature range. 
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Using this representation, the lattice specific heat is 



FIG. 10: Measured phonon density of states for T^MogSee, 
In2Mo6See and Rb^MoeSeg (histograms) plotted simultane- 
ously with PDoS data from neutron scattering (crosses). 

The fitted PDoS are shown for each compound in 
Fig. [lO] and the results obtained compare very favourably 
with previous neutron scattering dataP^ All three com- 
pounds clearly display the two key features already iden- 
tified in the C/T plots: a strong narrow peak below ~ 
10 meV and a broader distribution of phonons from ~ 
10-30 meV. The low-energy peak corresponds to an op- 
tical mode - the "guest" ion phonon - which is mediated 
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by the M ions vibrating in the tubes formed between the 
MogSeg chains. The broader maximum at higher energy 
is mainly due to acoustic intra-chain phonons, including 
the two Peierls modes discussed earlier. 

Comparing the PDoS of each compound measured, the 
two superconductors are very similar with a slight spec- 
tral weight shift to higher energy in the optical cation 
mode for In 2 MogSe6. This can be explained by consider- 
ing the smaller atomic mass of In compared to Tl. How- 
ever, the intensity of the optical mode in Rb 2 MoeSe6 is 
significantly reduced and its frequency has been shifted 
to ^8 meV, compared with ~4-6 meV in Tl 2 Mo6Se 6 
and In 2 Mo6Se 6 . Neutron scattering data also showed ex- 
tensive hybridization of the cation mode with the chain 
modes in Rb 2 MogSe6, thus inducing a deformation in the 
PDoS from 10-30 meV which can also be seen in our data. 



D. Normal-State Resistivity and Electron-Phonon 
Coupling 

The resistivity curves from 0.35-300 K for each com- 
pound are shown in Fig. [Tl] T^MogSeg and In 2 Mo6Se 6 
both show linear metallic behaviour in the normal state. 
There is no evidence for any negative curvature result- 
ing from strong correlation effects or the existence of a 
Luttinger Liquid ground state at high temperature.^ A 
saturation in the resistivity is observed for T < 15 K in 
Tl 2 Mo 6 Se 6 and In 2 Mo6Se6, with residual resistivity ra- 
tios of 10.1 and 8.8 respectively. Conversely, Rb 2 MoeSe 6 
undergoes a broad metal-insulator transition with a mini- 
mum at T c — 170 K and an activation energy Ea = 173K. 
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FIG. 11: Resistivity in M 2 Mo 6 Se 6 normalised to 300K. Inset: 
low-temperature resistivity in T^MogSee and h^MoeSee nor- 
malised to T 3 after residual subtraction, highlighting the large 
contribution from low-energy phonons in T^MoeSee- 

Using the same phonon frequency bins as those used to 
calculate the PDoS, we may now evaluate the electron- 
phonon coupling from the normal-state resistivity data 
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FIG. 12: Resistivity (with residual subtracted) normalised by 
T and fitted with the same Einstein mode energies used for 
the PDoS determination from our specific heat data. 



for Tl 2 Mo6Se6 and In 2 Mo6Se6 (it is not possible to per- 
form this analysis on Rb 2 MogSe6, due to the metal- 
insulator transition). This procedure has been success- 
fully used to obtain the electron-phonon coupling in the 
superconducting borides ZrBi 2 pD YB^H anc j j jU q 1 ^2\ 
as well as the clathrate superconductors Ba§Si46 and 
Ba 2 4Siioo, a more detailed account of the method can 
be found in the references. Our departure point is the 
generalised Bloch-Gruneisen formula!^ 



,„x , s 47TTO* 

Pbg(T) = p(0) + T 



Otj.-.fu} 



l) 2 dw 



(9) 



where x = u/T and al r F u is the electron-phonon 
"transport coupling function" which can be decomposed 
into Einstein modes to give 



,-FL = -z X! ^tr.k^kS(LO - LO k ) 



(10) 



Substituting this back into equation [9] yields the dis- 
crete version of the Bloch-Griineisen equation: 
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FIG. 13: Phonon density of states (histogram) and normalised 
electron-phonon transport coupling function at r F u (line) for 
Tl 2 Mo 6 Se 6 and In 2 Mo 6 Se 6 . 



p BG (T)=p(0) 



2tt 



Ea 



tr,k u k 



XkZ 



(e x * - l) 2 



(11) 



where the mode weighting parameters are the di- 
mensionless constants \t r ,k- The residual resistiv- 
ity ,o(0) is determined separately from the raw data, 
equalling 39.5 pfl cm for Tl 2 Mo6Se6 and 37.2 fj£l cm 
for In 2 Mo 6 Se 6 . 

Our fits are shown in Fig. [12] and display a clear dif- 
ference between Tl 2 Mo6Se6 and In 2 Mo6Se6- The initial 
slope of (R— Rq)/T is much steeper in Tl 2 Mo6Se6, due 
to a large contribution from a mode with energy 59 K. In 
contrast, the lowest energy mode contributing to the re- 
sistivity in In2MoeSe6 is centred at 103 K, with the first 
significant contribution only arriving at 180 K. In the 
absence of any data in the literature for the carrier den- 
sity n, we decompose the unscreened plasma frequency 
ne 2 /eom* and express our fitted values for a 2 
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normalised by ne 2 /m*. These are displayed in Fig 
superimposed on the PDoS. 

The principal electron-phonon coupling for each com- 
pound occurs in the 10-18 meV frequency window. How- 
ever, Tl 2 Mo6Se6 also exhibits a major coupling in the 
3.5-6 meV region, in direct contrast with In 2 Mo 6 Se 6 , 
which shows no coupling at all below 6 meV. cq r F u is 
intimately related to a 2 i^, the electron-phonon coupling 
function governing superconductivity, 35 ' 1 implying an ad- 
ditional electronic coupling to the low-energy guest ion 
phonon in Tl 2 MogSe6- 



IV. DISCUSSION 
A. Superconducting Transitions 

T^MogSeg and In 2 Mo6Se6 both display an anoma- 
lous broadening of their resistive superconducting transi- 
tions, whose amplitude and variation under applied field 
is not compatible with conventional 3D thermal fluctua- 
tion models. This suggests that the quasi-lD nature of 
these compounds has a significant influence on the size of 
the critical region around T c . In 3D systems, the broad- 
ening of a superconducting transition under applied mag- 
netic field is due to a finite size effect P3 the vortex- vortex 
separation limits the divergence of the correlation length 
at T c , hence reducing the coherence volume and increas- 
ing the importance of thermal fluctuations. In contrast, 
a perfect ID system cannot undergo a phase transition 
due to insufficient degrees of freedom. Tl 2 Mo 6 Se 6 and 
In 2 Mo 6 Se 6 lie in the crossover regime between these two 
extremes. 

Developing a theoretical model for this transition re- 
gion is a difficult task; however, we may consider the 
Peierls transition as an analogous crossover from a quasi- 
lD system to a quasi-3D ordered state. Theory predicts 
a suppression of T c by a factor of up to 4, together with 
a light smearing of the transition 3 -^ and certain Peierls 
systems indeed exhibit significantly broadened "jumps" 
in their resistivity as a result of quasi-lD fluctuations.^ 
It therefore seems reasonable to attribute the broaden- 
ing seen in Tl 2 Mo6Se6 and In 2 Mo6Se6 to the extreme 
low-dimensional nature of the compound. 

Explaining why the effect is so much more notice- 
able in Tl 2 Mo6Se6 than In 2 Mo6Se6, particularly in the 
specific heat jump, is rather harder especially given 
that In 2 MogSe6 is more anisotropic. The calculated 3D 
Ginzburg number G^d for Tl 2 Mo 6 Se6 is only 10% larger 
than that of In 2 Mo6Se6 and, in any case, G%dT c is sev- 
eral orders of magnitude too small to explain the ob- 
served broadening. Calculating the ID Ginzburg param- 
eters Gid from the Mishonov model yields more realistic 
transition width amplitudes, although the fact that the 
measured transition width for Tl 2 MogSe6 is larger than 
that of In 2 Mo 6 Se6 implies that this sample was less in- 
trinsically homogeneous. 

The evolution of Gid with applied magnetic field has 
not yet been calculated for a superconducting nanowire 
and we were hence unable to track the broadening of 
the resistive transition using a low-dimensional model. 
However, we expect that a field-induced finite size effect 
similar to that seen in 3D systems should control the 
transition widths. The crucial factor here is the differ- 
ence in coherence volumes £j_£\\ between the two super- 
conductors: 5290 nm 3 for Tl 2 Mo 6 Se 6 and 11400 nm 3 for 
In 2 Mo6Se6. It is well known that low-dimensional fluc- 
tuations play an increasingly important role upon the 
reduction of coherence length in a material; the smaller 
coherence volume in Tl 2 MoeSe6 must therefore outweigh 
its lower anisotropy relative to In 2 Mo6Se6. 
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It is clear that our understanding of low-dimensional 
fluctuations at a superconducting transition would 
greatly benefit from a detailed theoretical analysis. In 
particular, the reproducible deformation of the specific 
heat jump in Tl 2 Mo6Se6 (which presumably results from 
a displacement of states from above T c to the back of 
the jump) is a remarkable phenomenon and merits fur- 
ther attention. The shape of the jump is reminiscent 
of that seen in superconducting carbon nanotube ma- 
trices,^ suggesting that such behaviour may be intrin- 
sic to weakly-coupled superconducting filamentary net- 
works. A Berczinski-Kosterlitz-Thouless transition is 
also thought to occur in such materials. 40 Within this 
model, the superconducting transition is broadened due 
to the appearance of a phase-incoherent intermediate 
state (consisting of Josephson-coupled quasi-lD super- 
conducting fibres) separating the globally-coherent su- 
perconducting ground state at low temperature from 
the metallic normal state above T c . Tl 2 Mo6Se6 and 
In2MogSe6 would be prime candidates to undergo a sim- 
ilar transition. 

We note that despite the strong evidence for quasi- 
lD fluctuations around the superconducting transitions 
in Tl 2 Mo 6 Se 6 and In 2 Mo 6 Se 6 (as well as the broad re- 
sistivity minimum corresponding to the metal-insulator 
transition in Rb 2 Mo6Seg), there is no indication of a 
high-temperature Luttinger Liquid ground state in bulk 
Af 2 Mo 6 Se 6 . Transport experiments on MogSeg nanowires 
have revealed Luttinger behaviour in their conductance, 
which vanishes as the wire diameter increases above sev- 
eral tens of nanometers!^ However, since we are mea- 
suring M 2 Mo 6 Se6 crystals with diameters of the order of 
several hundred microns (and hence an increase in the 
number of conduction channels by a factor of 10 8 ), the 
absence of Luttinger effects is to be expected. 



B. Electron-phonon coupling 

Our BCS s-wave fits of the specific heat below T c dis- 
play conventional weak coupling (2Ao/fcsT c = 3.4) for 
In 2 Mo6Se6 and extremely strong coupling (2Ao/fcsT c > 
5) for Tl 2 Mo 6 Se 6 . In fact, Tl 2 Mo 6 Se 6 may well 
have usurped the throne of the /3-pyrochlore KOs 2 06 
(2Ao/fcsT c > 5) as the strongest-coupling phonon- 
mediated superconductor currently known. It should be 
noted that abnormally strong coupling (ranging up to 
2Ao/fcsT c ~ 25) is a characteristic of seve ral qu asi-lD 
CDW systems such as NbSe 3 and (TaSe 4 ) 2 l!^ISl due to 
the transition temperature being suppressed below its 
mean-field value. However, given that In 2 MogSe6 is more 
anisotropic than Tl 2 MogSe6, we do not believe that the 
strong coupling in Tl 2 Mo6Se6 originates from its low di- 
mensionality. 

Deconvolving the normal-state resistivity shows that 
the predominant common electron-phonon coupling for 
Tl 2 MogSe6 and In 2 Mo 6 Se6 lies in the 10-18 meV range, 
implying that superconductivity is mediated by the intra- 



chain modes which (according to neutron diffraction ex- 



periments^ 1 ) range from 12-40 meV and peak strongly at 
17 meV. This interpretation is supported by a tentative 
report of superconductivity under pressure in MoeSeg. 
The additional coupling to the low-energy optical mode 
in Tl 2 MoeSe6 moves this superconductor into the ex- 
treme strong-coupling regime, increases T c by nearly 2 K 
and reduces the coherence volume £^£n , rendering the su- 
perconducting transition more susceptible to broadening 
through quasi-lD fluctuations. 

Our observation immediately begs the question why 
In 2 MogSe6 does not enjoy a similar coupling to its op- 
tical In + mode. There are two reasons for this: firstly, 
consider the variation of the hexagonal lattice parameter 
a and the Pauling radii R p of the Tl + , In + and Rb + 
monovale nt ca tions. Values for a measured by X-ray 
diffractiorPE^in Tl 2 Mo 6 Se 6 , In 2 Mo 6 Se 6 and Rb 2 Mo 6 Se 6 



are given in Table IV together with standard R p values 
from the literature. Each cation is located at a/V3 from 
3 equidistant MogSeg chains and, since the chain radius 
is invariant with respect to the cation, the ratio \ZR p /a 
gives a good measure of the freedom of the respective 
cations to vibrate in their inter-chain tunnels. 



TABLE IV: 

M 2 Mo 6 Se 6 


Lattice parameters 


a and cation 


radii R p in 




a (A) 


flp(A) 


VSRp/a 


Tl 2 Mo 6 Se 6 
In 2 Mo 6 Se 6 
Rb 2 Mo 6 Se 6 


8.94 
8.85 
9.26 


1.15 
1.04 
1.48 


0.223 
0.204 
0.277 



It can clearly be seen that the In + ion is less geomet- 
rically constrained than Tl + and that Rb + is at consid- 
erably less liberty to vibrate than either of its "super- 
conducting" counterparts. This is evident in the neu- 
tron scattering data: as pointed out by Brusetti et at, 
Rb 2 Mo6Se6 displays significant hybridisation of the low- 
energy Einstein phonon, with the higher-energy inter- 
nal chain modes corresponding to a ^40% increase in 
the M ion force constants. Upon closer examination, 
the neutron-imaged PDoS of In 2 Mo6Se6 has a slightly 
deeper trough at ~11 meV than Tl 2 MoeSe6, implying 
marginally less phonon hybridisation, which is consis- 
tent with our estimate above. We therefore believe that 
the interchain tunnel diameter in In 2 MogSe6 is simply 
too large relative to the In + ion to allow its low-energy 
phonon to effectively couple to the Mo d electrons at the 
Fermi level in the chains. 

Secondly, the intrinsic electron-phonon coupling 
strength A is proportional to_l/oj 2 , where u is the char- 
acteristic phonon frequency.^ Due to its smaller mass, 
the In + mode is shifted to higher energy as can be seen 
both in our PDo S his tograms and the neutron data from 
Brusetti in Fig. 10 Using the Tl + = 5.2 meV and 
In + = 6.3 meV mode energies from Brusetti et al™, we 
calculate uj^/uj^ = 1.47; i.e. the coupling to the Tl + 
mode should be nearly 50% stronger than that to the 
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In + mode. To make a very c rud e comparison, we sum 
our measured a^F^ from Fig. 13 in the relevant energy 
range 5.1 - 8.9 meV, obtaining S a 2 F u (Tl)/T, a 2 F u (In) 
= 2.04. This suggests that the frequency-dependent vari- 
ation in coupling strength and the geometric constraints 
on the guest ion mode have a roughly equal importance 
in determining the coupling in M 2 Mo6Se6- 

Naively, we might expect an enormous electron-phonon 
coupling and ultra-strongly-coupled superconductivity in 
Rb 2 MogSe6 due to its narrow effective tunnel diameter. 
However, our electronic structure calculations indicate 
that it is not the geometric constraints on the M ion 
in M 2 Mo6Se6 which determine its anisotropy, but rather 
the degree of warping in its Fermi sheets. This is con- 
trolled by the in-plane bandwidth w, as detailed in Ta- 
ble |nj The electropositivity of the M ion - its willingness 
to donate electrons - is a useful quantity for character- 
ising the behaviour of Af 2 MoeSe6, since it is inversely 
proportional to w. Group IA metals are much more elec- 
tropositive than Group III and hence Rb 2 MoeSe6 falls 
victim to a high temperature insulating instability rather 
than becoming superconducting at low temperature. (It 
should nevertheless be noted that high-pressure measure- 
ments have succeeded in partially suppressing the metal- 
insulator transition in Rb 2 MogSe6 and simultaneously in- 
ducing superconductivity with a maximum T c ~ 4 K. 4 ) 
In a similar fashion, In is more electropositive than Tl, 
thus explaining the increased anisotropy seen in the su- 
perconducting state of In2MogSc6 compared to that of 
Tl 2 Mo 6 Se 6 . 



C. Metal-Insulator Transitions in A^MoeSee 

Our LDA calculations have indicated two possi- 
ble mechanisms for the metal-insulator transition in 
Af 2 Mo6Sc6: a Peierls transition resulting in the forma- 
tion of a dynamical CDW, or a SDW driven by strong 
Coulomb repulsion. Let us initially consider the SDW 
scenario. 

SDWs and their associated low-temperature antiferro- 
magnetic order are typically imaged using neutron scat- 
tering techniques. Unfortunately, no such data exist in 
the literature for any of the M 2 Mo6Se6 family. How- 
ever, a SDW may also be detected by the temperature 
dependence of its magnetic susceptibility x(T). In a 
SDW, x(T) generally exhibits paramagnetic behaviour 
above the transition, but then displays a characteristic 
maximum at the critical temperature Tmi before falling, 
hence signalling the onset of antiferromagnetism. Any 
SDW transition should therefore be clearly visible in our 
ac susceptibility data (Fig. [5]). Examining this closely, no 
evidence can be seen for any departure from temperature- 
invariant weakly-diamagnetic behaviour above the noise 
threshold of our susceptometer (10 -8 emu). In partic- 
ular, the ac susceptibility of Rb 2 Mo 6 Se 6 does not dis- 
play any anomaly as it passes through Tmi- We there- 
fore find no experimental evidence for SDW formation in 



M 2 Mo 6 Se 6 . 

Note that our calculated magnetic coupling constants 
A {/ are greater than the critical couplings A c for all 
M 2 Mo 6 Se6, not just Rb 2 MoeSe6. This implies that any 
SDW transition should be present in all three com- 
pounds. However, no Tmi may be identified in our 
resistivity data for Tl 2 Mo 6 Seg and In 2 Mo6Seg. Fur- 
thermore, it is not realistic to suggest that quasi- 
1D fluctuations suppress SDW formation in Tl 2 MoeSc6 
and In 2 MogSe6 but not Rb 2 MogSe6: from our elec- 
tronic structure calculations, Rb 2 Mo 6 Se 6 is much more 
anisotropic than Tl 2 Mo6Seg and In 2 MoeSe6 and should 
hence be more susceptible to such fluctuations. A far 
more probable scenario is that low-dimensional fluctua- 
tions prevent the SDW transition from taking place in all 
three compounds, with the metal-insulator transition in 
Rb 2 MogSe6 driven by a separate mechanism. Our LDA 
calculations for the electron-phonon coupling point to- 
wards this being a dynamical CDW and we will continue 
our discussion from this perspective. 

The transition seen in the electrical resistivity of 
Rb 2 Mo 6 Se6 is very broad and continuous compared with 
the distinct ste ps see n in other CDW materials such as 
NbSe 3 or TaS 2 P^Hl Such "blurring" can be explained 
by a combination of the influence of low-dimensional 
phase fluctuations and the dynamical nature of an un- 
derlying CDW whose order parameter may vary in both 
time and space. TarascorP and Hop^I have also per- 
formed transport measurements on Rb 2 MogSe6, obtain- 
ing similar broad transitions with Tmi ~ 125 K, 104 K 
and Ea = 87 K, 145 K respectively, which are rather 
smaller than our values Tmi ~ 170 K and Ea = 
173 K. We attribute the differences between Tmi to vari- 
ations in sample anisotropy, resulting from slight dif- 
ferences in the stoichiometries. It is well-known that 
Tl 2 Mo6Se6 is not a perfectly stoichiometric compound,^ 
with the highest Tmi occurring in Tli.gsMogSeg. Sim- 
ilar non-stoichiometric behaviour occurs in the rest of 
the M 2 MogSe6 family, for example producing occasi onal 
batches of In 2 MoeSe6 which are non-superconducting P32 

No clear evidence has been seen for the metal-insulator 
transition by any other experimental technique, although 
we observe a small anomaly in the heat capacity of 
Rb 2 Mo 6 Se 6 between 170 K and 185 K (see Fig. [8), just 
above Tmi- However, the height of the peak scarcely 
exceeds the noise threshold of our data and higher res- 
olution measurements will be required to verify and 
quantify this feature. Given the limited evidence for a 
thermodynamic phase transition at Tmi, it is tempting 
to hold intra-chain defects responsible for the insulat- 
ing behaviour observed at low temperature, relying on 
thermally-activated interchain hopping to enable metallic 
conductivity at high temperature. However, this simple 
model cannot be applicable to A/ 2 Mo 6 Se 6 since it is un- 
able to explain the broadband conductance noise indica- 
tive of sliding density wave motion!^- Furthermore, with- 
out a phase transition to deplete the DoS at the Fermi 
level Rb 2 MogSe6 would become superconducting at low 



14 



temperature (unless the defect density is sufficiently high 
to enable weak localization). No trace of superconduc- 
tivity in Rb 2 Mo6Se6 has been seen by any experimental 
technique at ambient pressure. 

An examination of the behavioural trends in the re- 
maining (Group IA) members of the A^MogSeg family 
lends further support to the CDW argument. Since the 
electronic anisotropy is proportional to the electropos- 
itivity of the M atom, we would expect to see metal- 
insulator transition temperatures fall as we move up 
Group IA (i.e. Cs—>Na). Although no systematic study 
of the family has been carried out by a single author, 
Tarascon et al™ found that Cs 2 Mo6Se6 has a higher Tmi 
than Rb 2 MogSe6, which is in line with this hypothe- 
sis. Hor et al™ have shown that hydrostatic pressure 
(leading to an enhanced interchain coupling and reduced 
anisotropy) suppresses Tmi to lower temperatures, im- 
plying that Tmi simply scales with the area of the Fermi 
surface available for nesting (as would be expected for a 
CDW). Nevertheless, a comprehensive study of all Group 
IA members of the M 2 MoeSe6 family will clearly be re- 
quired to validate our proposed dynamical CDW sce- 
nario. It is interesting to note that scanning tunnelling 
spectroscopy on small bundles of MogSeg chains has re- 
vealed metallic behaviour down to 5 KpD however the 
reduced electronic filling in MogScg results in a different 
electronic structure at the Fermi level with up to three 
bands occupied. : A simple density wave model is hence 
unlikely to be applicable in this compound. 

In addition to the obvious metal-insulator transition 
in Rb 2 Mo6Se6, it is important to address the possibility 
of superconductivity in Tl 2 Mo6Se6 arising from a hidden 
CDW ground state, as suggested by the Hall effect results 
from Brusetti et al.)^ Our results confirm that there is 
no discontinuity in the resistivity or jump in the spe- 
cific heat of either Tl 2 Mo 6 Se 6 or In 2 Mo6Se6 to support 
this conjecture. Local stoichiometric heterogeneity in 
the samples measured could produce a continuous series 
of local transitions, further blurred by quasi-lD fluctua- 
tions. However, this would produce a positive curvature 
in p(T) for T < 80 K, which is not seen. Alternatively, if 
we assume that the inter-chain coupling in T^MogSeg is 
sufficiently strong for such low-dimensional fluctuations 
to be negligible, a theoretical modeP^ may be invoked 
which defines two separate transition temperatures: T cu 
at which the CDW distortion occurs and T c i at which 
an energy gap opens over the entire Fermi surface and a 
metal-insulator transition occurs. Despite the presence 
of a structural modulation, metallic behaviour persists 
in the intermediate temperature range T c i < T < T cu 
whose width is acutely dependent on the anisotropy. Us- 
ing this model and the Brusetti Hall coefficient data, we 
could identify T - 80 K as T cu in Tl 2 Mo 6 Se 6 and as- 
sume that the superconducting transition takes place at 
T c > Td. However, we cannot justify disregarding quasi- 
lD fluctuations (as required by the model) when they are 
manifested so clearly in the deformation of the specific 
heat jump below T c in Tl 2 Mo6Se6- 



A more attractive explanation envisages a par- 
tial CDW gradually depleting the Fermi surface in 
Tl 2 Mo 6 Se 6 as the temperature decreases. This model 
agrees perfectly with the observed crossover from 
electron-like to hole-like carriers, yet a corresponding in- 
crease in carrier mobility below 80 K is still required 
to explain the linear metallic resistivity. It is not clear 
what type of physical mechanism could be responsible for 
such a rise in mobility and hence the reported Hall co- 
efficient data in Tl 2 Mo6Se6 remains mysterious. In con- 
trast, the situation in Rb 2 Mo 6 Se 6 is less complex: the 
observed transition to insulating behaviour implies that 
a gap opens across the entire Fermi surface. However, 
with the present data set it is unfortunately impossible to 
judge whether a similar "partial CDW" might be present 
for T > 170 K in this material. 

Recent worlP^ on (TaSe4) 2 I has suggested the exis- 
tence of so-called unconventional CDWs (UCDWs) in 
quasi- ID systems, with a fc-dependent gap which aver- 
ages to zero at the Fermi surface. This results in the for- 
mation of a pseudogap prior to a density wave transition 
and is supported by photocmission data in K0.3M0O3 and 
(TaSe 4 ) 2 lP 1 Within the pseudogap phase, the UCDW 
competes with low-dimensional quantum fluctuations as 
described by the Lee-Rice- Anderson modelP^ these act 
to suppress the Peierls transition temperature Tp well 
below its mean-field value Tmf, with the CDW en- 
ergy gap (the order parameter) fluctuating in time and 
space between these temperatures. It seems likely that 
Af 2 MogSe6 displays a non-trivial combination of low- 
dimensional fluctuation effects and a dynamical CDW 
with possible momentum dependence. Scanning tun- 
nelling spectroscopy (STS) or angle-resolved photoemis- 
sion spectroscopy (ARPES) would be ideal tools to verify 
the existence of an UCDW, fluctuating order parameter 
or dynamical CDW since although a pseudogap should 
open in the density of states, no static modulation of the 
structure or charge density is expected. 



V. CONCLUSIONS 

We have calculated the electronic structures and stud- 
ied the resistivity and specific heat in both superconduct- 
ing and normal states for the quasi-lD M 2 Mo6Seg family 
of cluster condensates (M=T1, In, Rb) . These materials 
lie on the border between superconducting and insulating 
instabilities. 

Superconductivity in Tl 2 Mo 6 Se 6 and In 2 Mo 6 Se6 is 
principally mediated by an internal phonon from the 
MogSeg chains with an energy in the 12-18 meV range. 
It is hoped that theoretical calculations will soon iden- 
tify the precise energy and nature of this phonon. 
Tl 2 Mo 6 Se 6 exhibits a further coupling to an optical mode 
at ^5 meV, which we identify with interchain vibrations 
of the Tl + ion. This extra coupling is not present in 
In 2 Mo 6 Seg due to the higher In + phonon energy and 
the reduced geometrical constraints on the In + ion be- 
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tween the MogSe6 chains. The low-temperature specific 
heat of Iri2Mo6See is well-fitted by a standard single-gap 
isotropic s-wave BCS model with 2A Q /k B T c — 3.5±0.1. 
However, the additional interaction with the low-energy 
mode in T^MogSeg pushes it into the strong coupling 
regime with 2Ao/fcsT c > 5. The specific heat jump at 
T c also exhibits a significant deformation with a shift in 
states to low temperature, which we attribute to strong 
low-dimensional fluctuations accentuated by the small 
coherence volume. STS or similar tunnelling experiments 
would provide conclusive proof of the gap symmetry, as 
well as displaying the coupling to the low-energy phonon. 
Our LDA calculations show that in all members of the 
M2M.OeSee family, a single ID helix band crosses the 
Fermi level. Its in-plane dispersion w is reduced by a 
factor of ten in Rh^MogSeg compared with Tl 2 Mo6Se 6 
and In 2 Mo 6 Se6, while the out-of-plane dispersion is prac- 
tically unchanged. Using an analytical model, we have 
shown that this reduction of w is sufficient to explain the 
trend from metallic conductivity followed by strongly- 



coupled superconductivity in Tl^MoeSee to a high- 
temperature metal-insulator transition in ftt^MogSeg. 
This insulating ground state is a consequence of either a 
dynamical Peierls (CDW) instability or a SDW transition 
driven by strong electronic correlations. The tempera- 
ture invariant magnetic susceptibility seen in all three 
compounds favours a CDW interpretation. A^MogSeg 
are therefore an ideal target for future ARPES experi- 
ments, since a dynamical CDW should generate a pseu- 
dogap at the Fermi level. 
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